Vibrational relaxation and dephasing of Rb 2 attached to helium nanodroplets 
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The vibrational wave-packet dynamics of diatomic rubidium molecules (Rb2) in triplet states 
formed on the surface of superfluid helium nanodroplets is investigated both experimentally and 
theoretically. Detailed comparison of experimental femtosecond pump-probe spectra with dissipative 
quantum dynamics simulations reveals that vibrational relaxation is the main source of dephasing. 
The rate constant for vibrational relaxation in the first excited triplet state 1 3 E^" is found to be 
constant 7 w 0.5 ns -1 for the lowest vibrational levels v < 15 and to increase sharply when exciting 
to higher energies. 



I. INTRODUCTION 

Helium nanodroplet isolation (HENDI) is a well- 
established technique for isolating molecules and form- 
ing clusters at low temperature (0.38 K) for spectroscopic 
studies Shifts and broadenings of spectral lines of 
molecules embedded in the helium droplets are small due 
to the weak dopant-host interactions as well as to the pe- 
culiar quantum properties of the superfluid helium nan- 
odroplets 0,01 ■ Nevertheless, the details of the solute- 
solvent interactions that are at the origin of the observed 
line shapes are currently being studied with various ap- 
proaches 041] • However, to date no time- resolved stud- 
ies aiming at resolving the details of the interaction of 
excited molecules with helium nanodroplets have been 
performed. 

Vibrational relaxation of molecules and molecular com- 
plexes embedded in helium nanodroplets has been stud- 
ied by the group of R. Miller using high-resolution in- 
frared spectroscopy and bolometric detection 0, lij-[TTj]. 
It was found that systems having a large energy gap be- 
tween the molecular vibration and the excitations of the 
helium {e.g. HF (v = 1)) couple very inefficiently to 
the helium environment, which leads to slow vibrational 
relaxation times t > 0.5 ms. The observed droplet-size 
dependent line shifts and broadenings point at coupling 
to surface excitations of the helium droplets (ripplons) 
being the main mechanism of relaxation pd| . Using mi- 
crowave spectroscopy, droplet-size dependent rotational 
relaxation times have been determined to be of order 1 - 
10ns 0. 

Recently, experiments have been performed that probe 
the dynamics of spin relaxation, dissociation, aggrega- 
tion, exciplex formation, and vibration of molecules at- 
tached to helium nanodroplets [l2| - [20j . These studies 
have in common that the helium nanodroplets act as a 
dissipative environment that decisively affects the out- 
come and the dynamics of the process of interest. There- 
fore, the dynamics of cooling, relaxation and dephasing 
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induced by the ultracold bath of helium atoms attracts 
an increasing amount of attention both from experiment 
and theory. 

The vibrational wave packet dynamics of alkali metal 
dimers attached to helium nanodroplets has been stud- 
ied using the femtosecond pump-probe technique in a 
series of experiments in our group 0, [l8l - l20j . Alkali 
metal atoms and molecules represent a particular class 
of dopant particles due to their extremely weak bind- 
ing to the surface of He droplets in bubble-like struc- 
tures [2ll - |24j . In particular, pump-probe measurements 
with K2 diatomic molecules in singlet states attached to 
helium droplets reveal a significant impact of the helium 
environment on the vibrational dynamics, suggesting the 
manifestation of a Landau critical velocity for the vibra- 
tional motion of K2 on the surface of superfluid helium 
nanodroplets fl8l l25|. 

In this work we present the detailed analysis of pump- 
probe measurements of the vibrational wave packet dy- 
namics of Rb2 molecules in triplet states attached to he- 
lium nanodroplets with regard to relaxation and dephas- 
ing induced by the helium environment. This system is 
particularly well-suited for a quantitative study of the 
molecule-helium droplet interaction due to the precise 
knowledge of the spectra and dynamics of gas-phase Rb2 
molecules and due to the weak molecule-helium interac- 
tions which allow for an accurate theoretical description. 

The long-lasting wave packet oscillations that we ob- 
serve up to delay times t > 1.5 ns suggested that the 
Rb2 molecules desorb off the helium droplets on a short 
time scale t < 10 ps and continue to vibrate freely in 
the gas- phase [2(J. This assumption was backed by the 
good agreement between the measured and theoretically 
predicted vibrational frequencies. Besides, the measure- 
ment of beam depletion using a separate detector was 
interpreted as clear evidence that excited Rb2 molecules 
desorb off the droplets on the time scale of the flight time 
from the laser interaction region to the detector (~ 1 ms) . 
Furthermore, earlier measurements with K2 dimers as 
well as theoretical simulations on K atoms attached to 
helium nanodroplets indicated desorption times in rang- 
ing 3 - 8ps and 10 - 30 ps, respectively [Tel [26j. 

A more detailed inspection of our data reveals, how- 
ever, that the Rb2 molecules are subject to continuous 
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vibrational relaxation due to the constant coupling to 
the bath of helium atoms on the time scale of the pump- 
probe measurements. The experimental signature of the 
coupling of vibrating Rb2 to the helium is the decreas- 
ing contrast of wave packet oscillation signals as well as 
changing amplitudes of individual Fourier frequency com- 
ponents due to the redistribution of populations of vibra- 
tional states. In particular, the pronounced dependence 
of the dephasing time on the quantum number v of ex- 
cited vibrational levels points at system-bath couplings 
being active. This observation is in line with earlier mea- 
surements of the fluorescence emissions of Na2 molecules 
in triplet states, which indicated vibrational relaxation 
in the excited electronic state to take place on the time 
scale of the life time of the excited state due to sponta- 
neous emission (~ 10 ns) [27|. Recently, it was observed 
that desorption upon electronic excitation may even be 
completely inhibited in the case of Rb atoms excited in 
a particular laser wave length range [28[ . 

Pioneering experiments on the vibrational dephasing 
and relaxation of molecules (I2) exposed to collisions 
with rare gas atoms at high density were performed by 
the Zewail group 12911 . m otivating theoretical studies by 
Engel, Meier et al. [30l [3l| . More recently, dephasing 
times as well as relaxation rates have been studied ex- 
tensively by the groups of Apkarian and Schwcntner by 
means of femtosecond spectroscopy of the vibrational dy- 
namics of halide molecules isolated in cryogenic rare-gas 
matrices (32T - I37) . Seminal work on time- resolved mea- 
surements of the dissipative fluid dynamics in bulk He- 
ll has been performed using femtosecond pump-probe 
spectroscopy of triplet He^ excimers created inside He- 
ll [38]. Due to the strong coupling of the highly excited 
He^ to the surrounding He which forms an extended bub- 
ble around Hcj the dynamics is fully damped after one 
period of motion. 

Theoretical studies on collisional quenching of rota- 
tions and vibrations of alkali dimcrs and other small 
molecules by helium atoms at low temperatures have re- 
cently been stimulated by the prospects of creating sam- 
ples of cold molecules using buffer-gas cooling as well as 
sympathetic cooling with ultracold atoms as a cooling 
agent [39T - l4l] |. In the system Li2+He, for instance, the 
quenching rate constants in the approximation of van- 
ishing temperature are predicted to increase by about 
one order of magnitude with increasing vibrational levels 

v = o- 10 m. 

In the case of weak couplings, which applies to our 
system, the concept of perturbations of the vibrational 
levels of the molecules by fluctuations in the bath modes 
is well-established. It leads to vibrational energy relax- 
ation along with the decay of vibrational coherences. 
There may well be additional pure dephasing mecha- 
nisms that are expected to vanish in the low temper- 
ature limit [iH, |43| . In the related energy gap picture 
the populations of individual vibrational energy levels v 
relax stepwise to the next lower vibrational energy lev- 
els v — 1 [44|. Here, the relaxation rate for an isolated 



vibrational level increases with the vibrational quantum 
number v [HI, EH EH . The evolution of coherences and 
populations of vibrational levels in the weak-coupling 
limit is often being modeled with the master equation 
description, obtained from the anharmonic molecular os- 
cillator coupled to a harmonic bath. Various coupling 
terms are used for describing different interaction mech- 
anisms [MEIE1-EII- 

Following established practice in the chemical literature, 
in this contribution we use "dephasing" to describe the 
general mechanism of loss of coherence between quantum 
states. Note that in other fields, this process is prefer- 
ably referred to as "decoherence" while "dephasing" is 
then used to describe that special occurrence of decoher- 
ence, where no dissipation is involved. 



II. VIBRATIONAL WAVE PACKET DYNAMICS 

The experimental arrangement used for recording fem- 
tosecond pump-probe photoionization transients is iden- 
tical to the one described previously [2(|. In short, a con- 
tinuous beam of helium nanodroplets of the size of about 
8000 He atoms is produced by expanding high-purity He 
gas out of a cold nozzle (Tss 17 K, diameter d = 5 /mi) at 
high pressure (p«50bar). The helium droplets are doped 
with two Rb atoms on average per droplet by passing 
through a pick-up cell that contains Rb vapor at a pres- 
sure pub ~ 2 x 10~ 4 mbar. Alkali atoms and molecules 
are peculiar dopants in that they reside in bubble-like 
structures on the surface of He droplets. Upon formation 
of a Rb2 diatomic molecule, the binding energy is dissi- 
pated by evaporation of helium atoms from the droplets 
and occasionally by desorption off the droplets of the 
newly formed molecule itself. This leads to an enrich- 
ment of droplet-bound Rb2 molecules in weakly bound 
triplet states. 

Further downstream, the doped He droplet beam in- 
tersects the laser beam inside the detection volume of a 
commercial quadrupole mass spectrometer. Due to the 
limited mass range, only bare Rb^ photoions are detected 
mass-selectively. The laser beam consists of pairs of 
identical pulses produced by a commercial mode-locked 
Ti: sapphire laser and a Mach-Zehnder interferometer to 
adjust the time delay between the pulses. The pulses 
have a duration of ~ 160 fs and a spectral bandwidth 
at half maximum of Aw; as ~ 80 cm^ 1 and peak pulse 
intensity ~ 4GW/cm 2 . 

A pronounced oscillatory photoionization signal is 
observed in the pump-probe transients for laser wave 
lengths in the range A = 960 nm - 1032 nm. A typical 
measured pump-probe transient recorded at A = 1006 nm 
is depicted in Fig. Q] (a). On the time scale of picosec- 
onds, the transient signal is modulated by wave packet 
(WP) oscillations with a period T% g rj 0.95 ps (see inset). 
In addition, this oscillation is amplitude- and frequency- 
modulated due to dispersion and subsequent revivals of 
the WP motion in an anharmonic potential. On the long 
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FIG. 1: Experimental (a) and simulated (b) pump-probe tran- 
sients of Rb2 formed on helium nanodroplets recorded at the 
laser wave length A = 1006 nm. 

time scale of hundreds of ps, the contrast of WP oscil- 
lations degrades monotonically and eventually vanishes 
at delay times > 1 ns. In this paper we argue that this 
slow decay of contrast is due to dephasing induced by the 
coupling of the Rb2 molecules to the helium droplets. 
Other sources of dephasing are conceivable: gas phase 
collision with evaporated gas and clusters of helium or 
dephasing due to the influence of rotations. However, an 
estimate of the Rb2-He gas phase collision rate gives a 
value far too small to account for the observed data. As 
for the influence of rotations, it is clear that an initial 
thermal population of rotational levels leads to a sim- 
ilar decay of signal contrast as observed in the experi- 
ment [491 ] - Again, in our case a detailed analysis shows 
that in order to account for the observed decay rates, un- 
physically large temperatures would have to be assumed. 
Moreover, the observed functional dependence of decay 
rates on the laser wavelength cannot easily be explained. 
Taking into account all these findings, our picture of vi- 
brational damping of the dimer through the interaction 
with the helium droplet allows for the most consistent 
explanation of all observed phenomena. 

A. Free gas phase dynamics 

Let us first briefly review the fundamental aspects of 
pump-probe spectroscopy of diatomic molecules isolated 
in the gas-phasc. A first pump pulse excites a coher- 
ent superposition of vibrational states (WP) in an ex- 
cited electronic state. After some time delay a second 
probe pulse projects the WP to a final ionic state which 
is detected as a function of time delay between the two 
pulses. Fig. [2] shows the relevant potential energy curves 
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FIG. 2: (Color online) Selected triplet potential-energy curves 
of neutral Rb2 and of the RbJ ionic ground state relevant 
to the present study. The arrows indicate the creation of a 
vibrational wave packet in the first-excited state followed by 
resonant two-photon ionization. 



of the Rb2 molecule in the triplet manifold. The straight 
arrows symbolize the pump and probe pulse excitation 
pathways. 

Due to the cold helium environment, at time t = only 
the vibrational ground state on the lowest triplet state 
a 3 E+ is occupied. Detailed experimental and theoret- 
ical analysis reveals that rotational degrees of freedom 
need not be taken into account (see also Sec. IIII[I . As 
outlined in Sec. Ill C[ we fully solve the time-dependent 
Schrodingcr equation involving all relevant potential en- 
ergy surfaces. Here, we want to point out that the main 
effect of the pump pulse is to generate a coherent su- 
perposition of vibrational eigenstates \v) on the excited 
electronic state surface (1) 3 S+. More specifically, the 

created WP can be written as \tp e (t = 0)) = J2v=o Cv \ v )> 
where c v denotes complex expansion coefficients and v = 
0,1,2,... is the vibrational quantum number. Coherent 
WPs will also be created in higher lying electronic states 
or, by resonant impulsive stimulated Raman scattering 
(RISRS) in the triplet ground state, with a significantly 
smaller amplitude, though [Ho) . Thus, all linear and non- 
linear processes are fully taken into account in our simu- 
lations. The respective populations of vibrational levels 
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|c„| 2 depend on the pump pulse parameters and on the 
Franck-Condon factors of the excitation transition. The 
pump pulse wave length determines the central vibra- 
tional level, while the pulse width, the pulse energy and 
the Franck-Condon factors determine the number and 
the relative populations of vibrational levels. The created 
WP in the (l) 3 E+-state propagates in the region between 
the classical inner and outer turning points. Direct inte- 
gration of the Schrodinger equation ihdt \ip e ) = H e \ip e ) 
yields 

N 

\i; e (t)} = J2^e-'^ 1 \v), (1) 

v=0 

where E v is the energy of the v-th vibrational level. As 
discussed below, the measured signal allows to extract in- 
formation about the density matrix of vibrational states 
p e {t). For the isolated dimer, p e (t) = \tp e (t)) (i(> e (t)\ de- 
scribes a pure state at all times. Its time evolution is 
given by 

Pe (i) = c v 4,e-' Ev ~° v ' * \v) (v'\ . (2) 

v ,v' 

The diagonal elements of the density matrix p vv = \c v \ 2 
represent the populations that are constant in time, while 
the off-diagonal elements p vv '(t) = (v\p e (t)\v') with v ^ 
v' oscillate with Bohr frequencies io vv i = (E v — E v i)/h and 
represent the coherences between the vibrational eigen- 
states \v), \v'). 

The probe pulse produces photoions through a reso- 
nant 2-photon-transition from the excited state to the 
ionic ground state 2 S+ of RbJ. Transitions preferably 
take place when the WP is located around a well-defined 
transition region, where the transition dipole matrix ele- 
ment is maximal (Franck-Condon window). Even though 
our simulations are numerically exact, it is instructive to 
consider the perturbative dependence of the ion signal on 
the WP density matrix [HHH^l , 

S(t) — ^ A vv > p vv > (t). (3) 

vv' 

The coefficients A vv i in Eq. @ contain products of tran- 
sition moments and field parameters and provide infor- 
mation about the vibrational populations in the final 
state. Through the dependence on the density matrix, 
the signal S is composed of beat frequencies ui vv i between 
all pairs of energy levels that contribute to the WP. The 
most prominent oscillation in the signal originates from 
components uj vv+ i and reflects the circulation of the WP 
on the potential energy surface. From the Fourier spectra 
of the signal information about higher-order frequency 
components uj vv+ ^ v with Aw > 1 can be extracted. In 
this way it is possible to gain information about the den- 
sity matrix from the measured ion signal. 

In a harmonic potential with energy levels E n = hu) e n, 
where n = 0, 1, 2, . . . denotes the number of cigenstate 



|n), Eq. ([T]) yields a periodic oscillation with classical 
period T c = h/ AE = 2ir/uj e , where AE denotes the con- 
stant energy spacing between adjacent levels n. The sig- 
nal S in Eq. ([3]) features a periodic oscillation with period 
T c . In the anharmonic Morse potential with energy spec- 
trum E v = fvjj e (y—x e v 2 ), however, initially well-localized 
WPs spread out due to dispersion on the characteris- 
tic time scale Tdi sp = Hlj / (uj e x e AE pump ) [37], HH, [56j . 
Here, x e is the anharmonicity constant, to denotes the 
central vibrational frequency of the WP and Ai5 pump 
is the spectral energy width of the pump laser pulse. 
Considering the parameters of our experiment we obtain 
Tdisp ~ 100 ps. Around t = Td; sp , all contributions on the 
right-hand side of Eq. ([3]) appear uncorrelated, which 
means that the oscillatory signal collapses There- 
fore, dispersion of the WP leads to a decay of the pump- 
probe signal contrast. Note, however, that due to disper- 
sion neither populations p vv nor the absolute values of 
coherences l/awl, v 7^ v', change in time. 

In the Morse potential, a revival of the initial WP 
takes place at certain times, i. e. the original phase cor- 
relation in the WP is restored and the WP partly or 
fully revives. At the full revival time, the signal ampli- 
tude ideally reaches its initial height, which underlines 
that coherence is preserved. Full revivals occur at times 
t = k x T rev /2, where T lev = 2ir/(uj e x e ), when all vi- 
brational cigenstates have accumulated a phase of 2itk 
with k = 1,2,3,.... At fractions of the revival time, 
t = p/qx T rGV where p/q is an irreducible fraction of inte- 
gers, the WP consists of a superposition of q copies of the 
original WP (fractional revivals) [H,[5^|. For instance, at 
half-period revivals (p/q = 1/2) the initial well-localized 
WP evolves into a highly quantum mechanical state that 
consists of two counter-propagating partial WPs that in- 
terfere with each other when colliding. In the electronic 
state of relevance for the present analysis, l 3 Sg of Rb 2 , 
which perfectly matches the shape of the Morse poten- 
tial in the accessible range of w-states, the first full revival 
time is T rcv /2 rj 160 ps [20|. 



B. Wave packet dynamics with dissipation 

The previous discussion was devoted to isolated vibrat- 
ing diatomic molecules in the gas-phase. Let us now 
consider Rb 2 molecules (M) coupled to the dissipative 
environment realized by helium nanodroplets (HND), to 
which the molecules are attached. 

A Rb 2 molecule attached to a HND is a closed but 
complicated system which can be described by the Hamil- 
tonian H = Hu + ^hnd + #m«hnd- Hm denotes the 
isolated molecule as discussed before, i?HND is the Hamil- 
tonian for the pure helium nanodroplct and -Hm-c+hnd 
contains the interaction between the two. As we are only 
interested in the dynamics of the molecule we call this 
our "system" and the helium nanodroplet our "bath". 
In our experiments the WP dynamics of the coupled Rb2 
molecule is mostly very similar to that in the gas phase, 
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which means that we see the same fast oscillations, WP 
dispersion with time constant ldi sp , and (fractional) re- 
vivals at times T rev . However, on the long time scale of 
the experiment (nanoseconds) the oscillatory signal ex- 
ponentially decays due to slow system dephasing. Thus, 
a description in terms of a weak system-environment cou- 
pling is justified. 

In the experiment, we observe a decay of the revival 
amplitudes with a rate 775. As we will show, this decay is 
related to the environment-induced dephasing of the WP 
due to dissipation. Dephasing can also be caused by a 
process which only affects the off-diagonal elements of the 
density matrix (no dissipation). This process is referred 
to as " pure dephasing" . For the well-known two- level 
system, the overall dephasing time constant T2 is related 
to relaxation (Ti) and pure dephasing without dissipa- 
tion (T») by 1/T 2 = 1/(2 Ti) + 1/T*. For multi-level 
oscillators, which we consider here, the relation between 
dephasing and dissipation is more subtle. No simple gen- 
eral expression connecting the corresponding time scales 
exists: depending on the shape of the WP, dephasing may 
take place on a much shorter time scale [60l - l63| . We re- 
call that it is important to distinguish between contrast 
decay due to dephasing - which is an irreversible pro- 
cess - and the reversible drop of the observed oscillation 
amplitude due to dispersion in an anharmonic potential. 

The dissipative vibrational dynamics is described us- 
ing the framework of Markovian master equations. At 
this stage, we do not aim at deriving such an equation 
from a microscopic Hamiltonian, which would require de- 
tailed knowledge of the helium "bath" and interaction 
Hamiltonians, -ffjjND and TJm^hnd, respectively. In- 
stead, we choose a well-established Markovian quantum 
optical master equation [(77l for a weakly coupled environ- 
ment. The density operator p(t) of the (reduced) system 
that describes dissipation in near-harmonic systems at 
zero temperature is given by 



dtp = tt 
in 



H M ,p\ + J2 ( U&i - \&\ L hP\ ) . (4) 



coupling to bath 

This equation is of Lindblad form [65[ . To describe dissi- 
pation, we use Li = y/^ca, where hi is the usual quantum 
mechanical ladder operator, defined through the harmon- 
ically approximated potential energy curve i. This Lind- 
blad operator induces vibrational dissipation on the time 
scale l/ji. More specifically, independently of the initial 
conditions, to good approximation, the mean energy of 
the excited WP decreases exponentially with a respective 
rate 7^. 

The relaxation rate constants 7, in Eq. (j4|) are taken 
as fit parameters to match the experimental data. The 
chosen dissipative Lindblad operator also affects the off- 
diagonal elements of the density matrix. The latter de- 
cay with time due to dephasing which implies a transi- 
tion from an initially pure state to a state mixture. For 
t — > 00, only the ground vibrational state is occupied, 
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FIG. 3: (Color online) First-order coherences (absolute val- 
ues of the off-diagonal elements of the density matrix, p v v +i), 
v = 0-6, for vibrational wave packet dynamics in the 1 3 S^"- 
state potential of Rb2 (b) in comparison with the harmoni- 
cally approximated potential (a). 



which is, of course, a pure state again. For the master 
equation Q it is well known that localized WPs are "ro- 
bust" in the sense that they suffer only little dephasing. 
In contrast, two such WPs separated by a (dimensionless) 
distance D in phase-space loose their coherence with an 
accelerated rate D 2 7 [6a469j |. 

This effect is visualized in Fig. El in which the time evo- 
lution of the first-order coherences |/3„„+i| of the dynam- 
ics in the Morse-type potential of the l 3 E+-state of Rb2 
(b) is compared with the dynamics in the harmonically 
approximated potential (a). The shown initial distribu- 
tion of coherences is obtained when exciting a WP formed 
of vibrations v ps 0-6 as observed in the experiment at 
a laser wave length A = 1025 nm. The matrix elements 
p vv > are computed by solving Eq. Q numerically and by 
projecting onto the cigenstatcs \v) at every time step. In 
the case of a harmonic potential (Fig. [3] (a)), beats be- 
tween low-lying levels n < 3 initially rise on the shown 
time scale, whereas the high-lying level beats monoton- 
ically fall to zero. This is due to the vibrational redis- 
tribution from the high-lying levels into the lower-lying 
ones, which evidently not only redistributes population 
but also transfers coherence from the high-lying to the 
lower-lying levels. 

In the case of the Morse- type potential (Fig.[3](b)), this 
tendency is much less pronounced. All the shown coher- 
ence terms drop down with slightly varying decay rates. 
The slower decay of the low-level beats (1, 2) and (2, 3) in 
contrast to the faster decay of (3, 4) and (4, 5) is partly 
reminiscent of the vibrational redistribution mentioned 
before. In addition, higher excited levels decay faster. 
Thus, for the Morse oscillator the overall dephasing ap- 
pears to be accelerated with respect to the harmonic os- 
cillator. Moreover, beats (2, 3) and (3, 4) in Fig.|3](b) are 



6 



periodically modulated with period T rev /2. In particu- 
lar, a slow loss of coherence or even a momentary increase 
of coherence is apparent at times close to the full vibra- 
tional recurrences when the WP is well-localized again. 
Around the half-period fractional revivals, when the WP 
splits into partial WPs that are maximally delocalized 
(large D), however, dephasing is fastest. 



C. Numerical simulation 

In order to reproduce the ion yield in the gas phase, we 
calculate the final state probability after the interaction 
with the laser field. For the isolated dimcr, we here follow 
the approach of [to| and fully numerically solve the time 
dependent Schrodinger equation 



:-ff M |*(t)) 



(5) 



for the full state vector ^ = -0 e , . . . ). The Hamilton 
operator Hyi now also contains the field interaction with 
the molecule. In particular, we take into account that 
the final state consists of the bound ion plus an ejected 
electron with energy E. Following the approach of (7l| . 
we use a discretization of the electronic continuum. We 
determine the final state probability \tpf(E, r)| 2 for dif- 
ferent pump-probe delays r and electronic energies E. 
Adding contributions with different energies E, we ob- 
tain a signal S(t), which is proportional to the gas 
phase ion yield. For the isolated molecule, Eq. ([5J also 
directly allows to obtain the density operator through 
p(t) = |$(f))(^ r (t)|. Potential energy surfaces and tran- 
sition dipole moments were provided by O. Dulieu [72j . 

For our phenomenological description of the helium in- 
fluence on the dimer dynamics, we switch to the density 
matrix description, Eq. Our aim is to ascribe cer- 
tain damping parameter values 7j(A) to the measured 
pump-probe signal at wave lengths A. From a numerical 
point of view, the evolution of the density matrix can 
become very costly, in particular, if one considers many 
potential energy surfaces and/or many vibrational states. 
We therefore return to an equation for the state vector, 
i. e. to a Schrodinger-type equation including relaxation 
(and thus, dephasing). The density matrix of the master 
equation Q is recovered from the stochastic Schrodinger 
equation on average, p(t) = \^ (t)) {9 (t)\ (73|. In practice, 
one has to determine many realizations of state vectors 
^SSE.ij which can then be used to extract the density 
matrix (Fig. [3|), coordinate, momentum, energy expecta- 
tion values , or the final state probability (Fig. H](b)) to 
compare with the experiment. 



III. PUMP-PROBE SPECTRA 

Upon laser excitation of Rb2 molecules formed on he- 
lium nanodroplets at wave lengths in the range A = 
960 nm - 1032 nm, coherent vibrational WPs are created 



in the first excited triplet state 1 3 S+ as well as in the 
lowest triplet state a 3 S+ by RISRS with varying rela- 
tive intensity. Around A = 1010 nm, the pump-probe 
signal as shown in Fig.[T](a) is dominated by WP motion 
in the l 3 E+-state. The amplitude modulation results 
from dispersion of the WPs and half-period as well as 
full recurrences arc observed with high contrast at revival 
times T rcv /4 ~ 80 ps and T rev /2 rj 160 ps, respectively. 
The nearly exponential decrease of the signal contrast is 
attributed to relaxation-induced dephasing and will be 
investigated in detail in the following. The simulated 
transient = 0.45 ns -1 ), depicted in Fig.Q](b), nicely 
reproduces both the vibrational recurrences as well as the 
overall damping due to vibrational dephasing. 

At short delay times t < 50 ps we observe slight devia- 
tions between the simulated and experimental transient 
signals, which we attribute to the influence of rotations, 
as seen previously for iodine molecules 7^. As men- 
tioned in the introduction, a dephasing influence of ro- 
tations on pump-probe spectra is conceivable j4^|, yet 
requires unphysically high temperatures > 10 K in our 
case. 

Gas phase simulations of the vibrational WP dynamics 
including free rotations in the Rb2 system feature notable 
rotational recurrences at half and full rotational periods 
T rot /2 w 575 ps and T rot ss 1150 ps, respectively. How- 
ever, no such signals are observed in the experimental 
data. Due to the surface position of the Rb2 molecules 
presumably with the molecular axis being oriented par- 
allel to the surface [75|], we expect the rotation to de- 
compose into weakly perturbed in-plane rotation and 
strongly hindered out-of-plane rotation which more likely 
resembles a pendular motion. The latter may efficiently 
couple to surface modes of the droplets causing fast relax- 
ation. Couplings between vibration, rotation, and libra- 
tion may therefore induce intricate relaxation dynamics, 
which, however, would require more expanded simula- 
tions that lie beyond the scope of the present work 

The experimental pump-probe signal of Fig. Q] (a) is 
analyzed by Fourier transforming the time trace inside a 
time window of a width of 5 ps and a Gaussian apodiza- 
tion function with full width at half maximum (FWHM) 
of 2.6 ps that slides across the data (wavelet or spectro- 
gram analysis). This type of analysis provides informa- 
tion about the frequency components that make up the 
WPs without losing the complete dynamical information. 
The result for A = 1006 nm is displayed in Fig. [4] (a). In 
this representation, the individual WP oscillations are 
no longer resolved, but the full, half-period and even 
one third-period revivals are clearly visible and can be 
attributed to frequency-beats between vibrational states 
separated by Av = 1, 2, and 3 vibrational quanta, re- 
spectively j59|. The fact that the WP recurrences are 
seen with such an extraordinarily high contrast even at 
long delay times is a consequence of the shape of the 
l 3 E+-potential that nearly perfectly matches that of the 
Morse potential in combination with weak system-bath 
couplings [20l | . 
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FIG. 5: (Color online) Exponential decay time constants ob- 
tained by fitting the maxima of the full (Aw = 1) as well as 
fractional (Av = 2, 3) revivals plotted against the average vi- 
brational quantum number of vibrational states initially pop- 
ulated by the pump pulse (v) . The filled symbols refer to the 
wave packet dynamics in the excited l 3 E^"-state, the open 
symbols refer to the a-state. 



FIG. 4: (Color online) Sliding window Fourier spectra (spec- 
trograms) of the pump-probe transients recorded at A = 
1006 nm. The width of the time window is 5ps in (a) and 
400 ps in (b). (c) displays the power spectrum of the integral 
experimental and theoretical transients. 



Fig. [4] (b) displays a magnified view of the spectro- 
gram of the same data when using a time window of 
width 400 ps and an apodization function with FWHM 
209 ps in the spectral range v = 34-35.5 cm -1 . In this 
representation of the data, the frequency resolution is 
comparable to the one obtained by transforming the in- 
tegral data set (Fig. 2] (c)) while still retaining the dy- 
namics on the long time-scale. The individual frequency 
components reflect beats between coherently excited ad- 
jacent vibrational states that are unequally spaced due 
to the anharmonicity of the potential. By comparing to 
the Fourier spectrum of the simulated data in Fig. @] (c) 
we conclude that the WPs excited at A = 1006 nm are 
composed of vibrational states v = 6-11 with relative 
amplitudes determined by the spectral intensity profile 
of the fs laser. 



A. Analysis of dephasing dynamics 

In a first attempt to analyze the loss of contrast of the 
coherent WP oscillation signal the envelopes of the beat 
signals corresponding to Aw = 1, 2, 3 are extracted from 
the spectrograms of the type shown in Fig. @] (a) at vari- 
ous laser wave lengths. Due to the low spectral resolution 
the few vibrational beats that are simultaneously excited 
according to the laser band width are not resolved. Thus, 
each envelope trace is composed of the sum of individual 
beats excited around the central laser wave length. The 
envelope traces are fitted to an exponential decay func- 
tion cx exp(— 7d^) to infer the characteristic decay rate 
7d- 

The fit results for 7d are depicted in Fig. [S] as a func- 
tion of the mean vibrational quantum number (v) deter- 
mined by the central laser wave length. The horizontal 
error bars reflect the width of the distribution of excited 
vibrational levels due to the laser band width, the vertical 
error bars depict the fit errors. Strikingly, the decay of 
contrast is strongly dependent on the level of vibrational 
excitation and features rapidly increasing decay rates 7d 
with increasing v. The solid lines in Fig. [5] (a) represent 
model curves obtained by fitting quadratic functions to 
the data. 

Similar behavior was observed in time-resolved cohcr- 
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ent anti-Stokes Raman-scattering measurements of the 
WP dynamics of molecular iodine I2 in the groundstate 
isolated in rare-gas cryo- matrices [33 - t34| . There, the 
transition from a linear w-dependence of 7 to a quadratic 
dependence with increasing temperature of the matrix 
was observed. Linear w-dependence at low temperatures 
was interpreted in terms of dephasing induced only by 
vibrational energy relaxation whereas at higher matrix 
temperatures pure elastic dephasing also contributed. 

Coherences between vibrational states spaced by Av = 
n > 1 are clearly subject to enhanced decay (Fig. [5] 
(a)). The experimental relative decay times amount to 

7d" =2 /7d , ' =1 ~ L3 and Td^/Td^ 1 ~ 2 as shown in 
Fig. O (b) . The scaling behavior of these multi-order de- 
cay rates 7d" = ™ with the order n of the beat has been 
discussed in the context of different mechanisms of pure 
dephasing [48[ . Depending on the collision model consid- 
ered in that study, a scaling behavior ranging from ze- 
roth to second order with n was expected. Our approach 
would require the inclusion of pure dephasing terms to 
account for Aw-dependent decay times whereas in the 
current model, where dephasing is solely induced by dis- 
sipation, the revival decay times are found to be inde- 
pendent of Av. 

The open symbols in Fig.[5]depict decay rate constants 
7g for the WP dynamics in the lowest triplet a 3 E+-state. 
While 7g of the first order coherence (Aw = 1) is similar 
to those of the excited l 3 S+-state for small (v), 7^ for 
the higher order coherences are significantly higher than 
for the l 3 E+-state dynamics at low (v). This result is 
reproduced by our assumption of vibrational relaxation 
and does not imply additional dephasing mechanisms. 
Note, that the initial distribution of w-state populations 
in the a-state is very different from that in the 1 3 S+- 
state. While in the l 3 S+-state several w-levels are popu- 
lated with similar intensities, in the a-state the popula- 
tion is peaked at v = and higher w-levcls are much less 
populated by RISRS. Therefore, the 7^ values to good 
approximation reflect dephasing rates between individual 
^-levels, since the beat signal of n-th order coherence is 
mainly composed of just one beat frequency. 

As we will see later, in the present case of Rb2 cou- 
pled to helium nanodroplcts, vibrational relaxation is 
likely to be the main source of dephasing. Yet, pure 
dephasing without population transfer does contribute 
to some extent. In a more complete description both 
dissipative as well as additional pure dephasing terms 
should be included to account for the observed v- and 
Aw-depcndcnccs. For the sake of restricting the model 
to the essential features of the problem, however, in the 
following discussion we focus on the model calculations 
that are based on vibrational relaxation. 



B. Numerical simulation 

In order to obtain a more quantitative description of 
the observed dynamics, the experimental data are mod- 



eled using the method outlined in Sec. Ill CI The only 
adjustable parameters entering the simulation are the en- 
ergy relaxation rate constants in the triplet ground and 
first excited states, j a and 7s s , respectively, as well as 
relaxation constants for the two probe states 3 3 n„ and 
4 3 E+. The resulting spectrograms of the best fits to the 
experimental data are displayed in Fig. [6] (bottom row) 
for the selected laser wave lengths A = 1025, 1006, and 
970 nm. For comparison, the top row depicts the exper- 
imental data and the middle row shows the simulation 
when relaxation is absent (70 = 7s g = 0). 

The transient recorded at A = 1006 nm (center col- 
umn in Fig. [6j , already shown in Fig. [H is dominated by 
the fundamental as well as by the first and second over- 
tone beats of the l 3 S+-state. The experimental data 
(Fig. [6] (d)) are very well reproduced by the numeri- 
cal simulation for a damping constant j-v, g = 0.45 ns _1 
(Fig. [5] (f)), whereas the agreement is clearly worse when 
no damping is assumed (Fig. [6] (e)). At laser wave 
lengths A = 1025 nm and A = 970 nm, WP oscillations in 
both ground a 3 E+ and excited states 1 3 S+ are present. 
At A = 1025 nm, the excited state-dynamics clearly 
fades away more slowly than at A = 1006 nm, which 
is in agreement with the simulated data when setting 
7 Sg =0.36ns~ 1 (Fig. EI (c)). In contrast to the l 3 S+-state 
WP-dynamics, the a-state beats feature less visible dis- 
persion and recurrences of the WP motion. This is due to 
the fact that predominantly the vibrational groundstate 
V = is populated by RISRS. Consequently, the funda- 
mental spectral component uj a ~ 13 cm -1 is mainly com- 
posed of the beat frequency (E v= \ — E v = )/(hc), with lit- 
tle contributions of (E v= 2 — E v —i)/(hc) and higher level 
beats [20I ] . Best agreement with the experimental data is 
obtained for j a = 3ns _1 (Fig. [6] (c)). When no damping 
is assumed, the simulation clearly severely deviates from 
the experimental data (Fig. [5] (b)). We attribute the ad- 
ditional spectral features to the dynamics in the higher- 
lying electronic state 3 3 II U which has a similarly shaped 
potential curve as the l 3 S+-state. In all simulations in- 
cluding vibrational relaxation, the damping constants of 
the 3 3 n„ and of the 4 3 S+-states are set to O.lps -1 to 
achieve fast damping of the corresponding WP dynamics, 
no direct WP signal related to these states is observed 
in the experiment. However, the inclusion of the two 
states is crucial in order to reproduce the experimentally 
observed WP signals in the a-state. The excited state dy- 
namics at A = 970 nm is only visible in the time range 
- 20 ps, after which the a-state dynamics prevails. This 
behavior is reasonably well reproduced by the simula- 
tion when assuming very fast relaxation (js g = 0.5 ps -1 , 
Fig- [HI (i))- ^ n contrast, the same simulation with 7s g = 
shows a dominant contribution of the excited l 3 E+-state 
dynamics (Fig. [6] (h)). At A = 980 nm (not shown in 
Fig. [S]), l 3 S+-state components are still visible during 
delay times - 100 ps, which implies fast relaxation at a 
rate jz g = 0.01 ps -1 . 

An even more detailed verification of the numerical 
model is achieved by comparing the experimental and 
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FIG. 6: (Color online) Comparison between experimental (top row) and theoretical (middle and bottom row) data in spec- 
trogram representation at selected laser wave lengths; the middle row shows the simulation of undamped vibration without 
coupling to the helium droplets; the bottom row shows the simulation including damping. 
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FIG. 7: (Color online) Time evolution of individual beats 
between adjacent vibrational levels (Av = 1) of the l 3 £jj~- 
state extracted from spectrogram analysis of the transient 
at A = 1025 nm with a 400ps-time window ((a) and (b)) in 
comparison with simulated density matrix elements |p„„+i| 
((c) and (d)). Plots (b) and (d) show the same data as (a) 
and (b) where the beats are normalized to the sum of all beats 
at each time step. 



theoretical data in the spectrogram representation using 
a long time window of width 400 ps, as shown in Fig. [4] 
(b). The high spectral resolution retained in this anal- 
ysis allows to compare the time evolution of individual 
beats between adjacent vibrational states. The ampli- 



tudes of individual frequency components are extracted 
from vertical cuts through the spectrograms at maximum 
positions and are plotted in Fig.[7Jfor A = 1025 nm. Pan- 
els Fig. [7] (a) and (b) represent the experimental data, 
where in (b) each amplitude component is normalized to 
the sum of all contributing beat amplitudes. Although 
all of the frequency components except the lowest one 
(v = 0, v = 1) decay in absolute amplitude (Fig. [7] (a)), 
the relative amplitudes only decrease in the case of the 
high-lying level beats (4, 5) and (5, 6), whereas the lower 
beats (3, 4) remain constant or even rise [(0, 1), (1, 2), 
and (2, 3)] in amplitude in proportion to the sum of all. 

The numerical simulations (Fig. [JJ (c) and (d)) show 
the evolution of the first-order coherences of the density 
matrix. The good agreement justifies the assumed model 
based on vibrational relaxation and highlights the pos- 
sibility of extracting information about the density ma- 
trix by appropriately analyzing the measured ion yields. 
The general decay and oscillatory behavior of individ- 
ual beats is well reproduced by the numerical simulation 
except for the (0, l)-beat which is extraordinarily promi- 
nent in the experimental data. A slight increase of the 
absolute beat amplitude of the (0, l)-component can only 
be explained by a redistribution of population of higher- 
lying vibrational levels into lower-lying ones. The weak 
periodic modulations of both experimental and theoret- 
ical curves are reminiscent of the revival structure that 
becomes more pronounced as the Fourier time window 
is reduced. The different decay rates for the individual 
vibrational beats have been discussed in terms of vibra- 
tional redistribution in the harmonic and anharmonic os- 
cillators, (Sec. HHland Fig.©. 

Since all of the beats are subject to dephasing the vi- 
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FIG. 8: (Color online) (a) Time evolution of the populations 
of individual vibrational states extracted from the simulation 
at A = 1025 nm. (b) Evolution of the vibrational level v pop- 
ulated on average and of the mean vibrational energy (E v ). 

brational redistribution is masked by an overall decay in 
the representation of absolute amplitudes in Fig. @] (a) 
and (c). However, from the simulation we can extract 
information about the evolution of populations of the in- 
dividual vibrational states. To this end, the diagonal el- 
ements of the density matrix are computed for each time 
step by projecting the wave function onto the vibrational 
eigenfunctions. The resulting populations of levels v = 
- 6 are depicted in Fig. |8] (a). Fig. [8] (b) shows the time 
evolution of the quantum number of the vibrational state 
that is populated on average as well as the corresponding 
average vibrational energy. Accordingly, at A = 1025 nm 
the vibrational populations relax down by about 1.8 vi- 
brational quanta during 1.5 ns. The corresponding vi- 
brational energy is reduced by about -Ediss = 73 cm -1 . 
Note that at shorter wave lengths the amount of de- 
posited vibrational energy into the droplets in this time 
interval is considerably larger, e.g. i?diss = 157 cm -1 at 
A = 1006 nm and -Ediss = 656 cm -1 at A = 970 nm. 

At such high rates of energy transfer to the helium 
droplets one has to consider the droplet response in terms 
of heating, superfluidity and cooling by evaporation of 
helium atoms. Dissipation of vibrational energy up to 
-Ediss = 656 cm -1 into the droplets leads to a significant 
rise in droplet temperature, which may locally exceed the 
transition temperature to the supcrfluid phase (2.17 K for 
bulk helium) of even the boiling point. Note, however, 
that cooling of the droplets due to evaporation of he- 
lium atoms may counteract the heating process. In the 
considered energy range, effective cooling is expected to 
set in on a time scale of ~ 100 ps [76|. Thus, slow en- 
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FIG. 9: (Color online) Dependence of the damping parameter 
obtained from fitting the simulation to the experimental data 
as a function of the average vibrational level populated at 
different laser wave lengths. Filled and open symbols refer to 
the 1 3 S^"- and the a-states, respectively. 



ergy transfer from the molecules to the droplets at excita- 
tions to low-lying ^-levels could be partly compensated 
by evaporation of helium atoms (Ediss ~ — 5cm _1 per 
evaporated atom) , whereas fast energy transfer may lead 
to effective heating, to subsequent local disequilibrium 
states and even to the breakdown of superfluidity. 

However, after the time of flight of the droplets to 
the beam depletion detector (~ 1 ms) the complete vi- 
brational excitation energy will be transferred to the 
droplets leading to massive evaporation of up to 150 he- 
lium atoms. Thus, the observed beam depletion signal 
could be due to the deviation of the evaporating droplets 
out of the beam axis instead of being the result of des- 
orption of Rb2 from the droplets, as assumed previously. 
This assumption is supported by the observed significant 
blue-shift of the beam depletion signal with respect to the 
pump-probe-photoionization spectrum (Fig. 11 in [2fj|). 
In the spectral range A > 1000 nm, in which a high pho- 
toionization yield is observed but only slow dissipation, 
beam depiction is low. 

Having established the observed decay of vibrational 
beats in terms of relaxation-induced dephasing using our 
numerical simulation, let us finally discuss the relaxation 
time-constants 7, obtained by fitting the model to the ex- 
perimental data. Fig. [9] shows j a and 7e 9 as a function 
of the average vibrational quantum number correspond- 
ing to WPs created at different laser wave lengths (top 
scale). Interestingly, in the range v = 2 - 14 the damp- 
ing parameter remains nearly constant, 7s 9 fts 0.5 ns , 
even though significantly varying decay rates 7d have 
been measured (Fig. [SJ. This discrepancy reflects the 
scaling behavior of dephasing times with v, as discussed 
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in Sec. Ill Bl The nearly u-independent values of js g are 
quite unexpected considering various model predictions 
of strongly w-dependent relaxation rates [H, 0, IZZt ■ At 
higher vibrational excitations v > 15, though, the exper- 
imental data can only be modeled when assuming drasti- 
cally increased values of 7s s . The ground state relaxation 
rate j a (open symbol in Fig. [5]) is found to be higher by 
about a factor 6 as compared to the excited state rate 

The nearly constant low values of in the range of 
small values of v may be related to significant mismatch 
between energy spacings (~ 35 cm -1 ) and excitation en- 
ergies of the helium bath modes. The ripplon energies 
are in the range 0.1 cm' 1 , phonon modes have energies 
~ 1 cm -1 , and the roton energy is about 10 cm -1 [76Ll78j |. 
Consequently, the ground state vibration with lower level 
spacing (~ 13 cm -1 ) more efficiently couples to the he- 
lium environment, which would explain the higher re- 
laxation rate j a . The sharp rise of 7s g at (v) > 15 
could be related to the breakdown of superfluidity or 
even to the effect of a liquid to gas phase transition at 
the Rb2-He interface induced by fast heating. New cou- 
pling channels, e. g. the excitation of collective modes of 
the helium droplets that may be related to their super- 
fluid character (i. e. rotons) may also be at the origin 
of increasing relaxation rates |2{|. At this stage, how- 
ever, this assumption seems unlikely, given the fact that 
the vibrational energy quanta (~ 35 cm ) largely exceed 
the elementary excitations of supcrfluid helium droplets. 
Possibly, intra-molccular couplings or more complex ex- 
citation pathways leading to the ionic continuum may 
also be involved. 



IV. CONCLUSION 

In conclusion, femtosecond pump-probe measurements 
of the vibrational wave packet dynamics of Rb2 molecules 
attached to helium nanodroplets are analyzed using dis- 
sipativc quantum simulations. In contrast to earlier in- 
terpretations, Rb2 excited to triplet states are found to 
remain attached to the helium droplets on the time scale 
of the pump-probe experiments and reveal slow damping 
of the vibrational wave packet signal due to the inter- 
action with the helium droplet environment. The weak 
system-bath coupling results in slow damping dynamics 



compared to the periods of vibration, which is proto- 
typical for the applied master equation. Thus, helium 
droplets provide a versatile test bed for studying relax- 
ation dynamics and cooling processes induced by a highly 
quantum environment. 

From the detailed comparison of the experimental data 
in the time- and frequency domains with model calcula- 
tions it is possible to deduce the evolution of the density 
matrix describing the vibrating Rb2 molecules. While 
rotational-vibrational coherences as well as pure dephas- 
ing of the vibrational wave packet dynamics may play a 
minor role, good agreement is achieved under the model 
assumption of vibrational relaxation-induced dephasing. 
We extract damping constants for the vibrational relax- 
ation in the lowest triplet state a 3 E+ and in the first 
excited state 1 3 S+, j a w 3ns _1 and ~ 0.5ns -1 , re- 
spectively. The pronounced dependence of "f^ g on the 
vibrational quantum number v may be related to the in- 
terplay of effective heating of the droplets by fast energy 
transfer from the molecules and cooling due to evapo- 
ration of helium atoms. High heating rates may induce 
phase transitions in the droplets that affect the dynamics 
of attached molecules. 

Further experiments as well as modeling of the re- 
sponse of the helium droplets to vibronic excitations of 
embedded atoms and molecules are needed. In particu- 
lar the dependence of dephasing and relaxation dynam- 
ics of the vibrational state quantum number may pro- 
vide detailed information about the solute-solvant cou- 
pling mechanisms and may contribute to interpreting the 
shapes of spectral lines. Moreover, the dynamics of the 
desorption process of alkali atoms and molecules off the 
droplet surface in dependence of the atomic species and 
the particular vibronic state excited is still mostly unre- 
solved. 
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